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Abstract 

Isospectral domains are non-isometric regions of space for which 
the spectra of the Laplace-Beltrami operator coincide. In the two- 
dimensional Euclidean space, instances of such domains have been 
given. It has been proved for these examples that the length spec- 
trum, that is the set of the lengths of all periodic trajectories, coin- 
cides as well. However there is no one-to-one correspondence between 
the diffractive trajectories. It will be shown here how the diffractive 
contributions to the Green functions match nevertheless in a "one-to- 
three" correspondence. 
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1 Introduction 



The quantum-mechanical problem of finding isospectral domains, that is two 
non-isometric regions for which the sets {En, n e N} of solutions of the 
stationary Schrodinger equation 

{A + E)^^0, (1) 

with ^ I boundary = 0, are identical, has been formulated in a synthetic way 
by Mark Kac in 1966 in his famous paper "Can one hear the shape of a 
drum" [1]. Negative answers to this problem have been given for domains 
on Riemannian manifolds [2, 3]. The answer for two-dimensional Euclidean 
domains was finally given by Gordon et al. [4], who provides explicitly a 
pair of simply connected non-isometric Euclidean isospectral domains. The 
two billiards considered in [4] are represented in Figure la. Using a paper- 
folding method. Chapman [5] showed that in this case isospectrality arises 
from the existence of a map between the two domains. This method allowed 
him to produce more examples of isospectral billiards (by billiard we mean 
a two-dimensional Euclidean connected compact domain): he showed that 
any triangle or even rectangle could replace the base right angled isosceles 
triangle used to build the billiards in la; one just has to glue together 7 copies 
of the chosen base shape obtained by symmetry with respect to its edges, 
in the same order as to build the biUiards of la. We will work here with 
the pair of billiards Bi and B2 of Figure lb, where each billiard is made of 
seven rectangles. Later on, Buser et al [6] produced more examples of planar 
isospectral domains made of more building blocks; all of them are based on 
the same principle of gluing together copies of a base triangle. 




Figure 1: a) Two isospectral billiards with a triangular base shape; b) The 
same with a rectangular base shape. 

A natural question arises when one considers isospectral billiards: is there 
any relation between their periodic orbits, and is there any relation between 
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their diffractive orbits? It is well-known that the quantum density of states 



(2) 



n 



can be expressed by means of the advanced Green function as 



d{E) = Im / daG{a,a), 



(3) 



where the integral is performed over the domain. In the case of two isospectral 
billiards, one might naturally expect the integrals of the Green functions of 
the two systems to be equal. Since the Green function can be expanded 
as a sum over all Feynman paths, the correspondence between the spectra 
should be associated to a correspondence between the periodic orbits and 
between the diffractive orbits. Moreover, there exists an exact expression 
for the Green function in a two-dimensional polygonal billiard: following 
Sommerfeld [7, 8], who provides the exact Green function for a half infinite 
straight mirror in 2 dimensions, Stovicek [9, 10] expressed the Green function 
for a collection of magnetic flux lines on a plane (the multi-fiux Aharonov- 
Bohm effect) as a sum over all possible scattering paths. This method has 
been generalized in [11] to provide the exact Green function for the scalar 
wave equation in a plane with any set of perfectly reflecting straight mirrors 
joined by diffractive corners. The Green function is given as a scattering 
series involving all classical trajectories and all scattering contributions; a 
semi-classical series expansion shows that one expects the Green functions of 
the two system to correspond order by order. 

The correspondence between the periodic orbits of two isospectral bil- 
liards (which contribute to the lowest order of the scattering series) has been 
discussed in [12] for isospectral domains on Riemannian manifolds. The 
Laplace spectrum versus the length spectrum is also discussed by Gordon in 
[4] for Euchdean isospectral billiards. A proof of the one-to-one correspon- 
dence between the length spectra (referred to as "iso-length spcctrality" ) in 
the case of the two celebrated isospectral billiards considered by Gordon is 
given in [13], based on simple mathematical tools. More generally, iso-length 
spectrality has been proved in [18] for all pairs of isospectral billiards hav- 
ing a "transplantation" property (which corresponds, roughly speaking, to 
the existence of a map between the two domains, such as the one that will 
be given in Section 2.2). However, what [13] underlines is that there is no 
one-to-one correspondence between the diffractive trajectories of the two bil- 
liards. More precisely, one can find straight lines between diffractive corners 
in one of the billiards that do not have any counterpart in the other. Still, 
the equality of the spectra, and hence the equahty of the trace of the Green 
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function, imposes some correspondence. In Section 2 we will study the map 
that exists between the two billiards and show how this bijection allows to 
prove a one-to-one correspondence of the periodic orbits. In Section 3 we 
establish a relation between the Green function of the two studied domains. 
Finally in Section 4 we answer the following question: what is the correspon- 
dence between the diffractive orbits. The conclusion briefly shows how the 
method presented here can be naturally extended to other pairs of isospectral 
biUiards. 

2 Isospectrality and periodic orbits 
2.1 The translation surfaces 

Instead of studying directly the biUiards Bi and B2 of Figure lb, we will 
consider the equivalent problem of studying the translation surfaces [14] as- 
sociated to these billiards. Bi and B2 arc polygons with angles (7r/2, n/2, 
7r/2, 7r/2, 7r/2, 7r/2, 7r/2, 37r/2, 37r/2, 2n). A construction due to Zemlyakov 
and Katok [15] shows that the translation surface associated to a generic 
rational polygonal billiard is obtained by unfolding the polygon with respect 
to each of its sides, which means gluing to the initial polygon its images by 
reflexion with respect to each of its sides and repeating the operation. If 
ttj = -KiTLi/ni are the angles of the polygon and N the least common multiple 
of the Hi, then 2N copies of the initial billiard are needed. In the case of the 
biUiards Bi and B2, since all the angles are multiples of 7r/2, only 4 copies 
are needed, and the translation surfaces M.i and M.2 obtained by this con- 
struction are represented in Figure 2. In these surfaces, all opposite sides are 
identified. Note that each structure has 4 singularities: 2 angles of measure 
Gtt (dot and circle), and 2 angles of measure 47r (cross and star). Moreover, 
each singularity with angle 2{k + 1)ti brings the contribution k to the quan- 
tity 2g — 2, where g is the genus of the surface [16]. Therefore the surfaces 
Ail and M.2 both are of genus 4. Any path drawn on Aiv corresponds to an 
unfolded path on By. 

Each surface is made of 7 tiles, which give to each translation surface the 
structure of a 7-fold torus cover. If the tiles are numbered as in Figure 2, 
the way these tiles are glued together to form a surface of genus 4 can be 
expressed, following [13], through three 7x7 matrices: for each structure 
AAu-i V = 1,2, we introduce the matrices R'^'^\ U^^'^ and U^'^'' such that 
Rf- = 1 if the right edge of tile i is glued to the left edge of tile j and 

otherwise, ?7 • = 1 if the left half of the upper edge of tile i is glued to the 
left half of the lower edge of tile j and otherwise, and of course uf^ = 1 if 
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Figure 2: 4 copies of the two isospcctral billiards Bi and B2 are glued together 
to make two isospectral translation surfaces A^i and AI2 niade from 7 tiles. 

the right half of the upper edge of tile i is glued to the right half of the lower 
edge of tile j and otherwise. Just looking at Figure 2, we can obtain these 
matrices easily. They are given in Appendix A. For instance, in Ali, the 
right neighbour of tile 5 is tile 1, the right neighbour of tile 3 is tile 7, etc... 
We can also define the matrices L^'^^ , D ('^^ and D ('^^ which indicate which tile 
is glued to the left, the bottom left or the bottom right of a given tile. These 
three matrices are nothing but the transposes of respectively R^'^\ U^'^^ and 

2.2 Mapping between translation surfaces 

The isospectrality of the two billiards arises from the existence of a mapping 
between A4i and A^2, provided by Gordon [4] and made explicit by Chapman 
[5]. The following section proves the isospectrality between M.i and A42- 

By convention, we will label any point a in or (which are made 
of 7 tiles ore tori) by its position a in the tile and the number i (1 < i < 7) 
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of the tile it is in; we will write alternatively a or (a,i). Let us define the 
"transplantation matrix" T as 
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(4) 



The isospectrality arises from the fact that for any given eigenstate 0„ of 
All, we can construct an eigenstate ipn in M. 2 defined by 

'4^n{a,i) ^ ^^Tij(pn{a,j), (5) 

where An is a normalization factor that we will discuss in Section 3. For 
instance, for a in tile 1 of M.2, we have 

1) = ^ i<Pn{a, 1) + 0n(a, 4) + (f)„ia, 7)) . (6) 

We call the tiles 1,4,7 in Aii the "pre-images" of tile 1 in A^2, and we say 
that 1 is "made of" tiles 1,4 and 7. The fact that the functions ipn are the 
eigenstates of comes from the following relations: 

jimT^TR^^\ U^^^T ^TIJ^^\ U^^^T ^TU^^\ (7) 

which can be verified directly on the matrices given in Appendix A. Each 
of these commutation relation has a natural interpretation. For instance the 
commutation relation for R means that for any pair {i-ij) of tiles, i on AI2 
and J on All, we have 

^2 = X] ' (8) 

k k' 

which means that if k is the tile on the right of tile i in Ai2 (i.e. uf^ = 1) 
and if we call ii,i2,i3 the pre-images of i (i.e. Ta^ — Ta^ — Ta.^ = 1), then 

by (8) 

+ + (9) 

Among the 7 possible values of j, the right-hand side of Equation (9) will 
be 1 if and only if j is on the right of ii, 12 or in Ali, and otherwise. 
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This implies that if k is on the right of i, the three pre-images ki, k2, k^ of k 
(verifying T^^ki = 1) are the three tiles on the right of the three pre-images 
of I. Concretely, since tile 1 in M.2 is made of tiles 1,4,7, then its right 
neighbour tile 5 is made of the right neighbours 2,4,6 of tiles 1,4,7. Since the 
commutation relation (7) is valid for all matrices A G {R,L, U , U , D , D} 
all properties of continuity between tiles in A^i are preserved by transplan- 
tation. Therefore all the functions ipn constructed by (5) are continuous on 
M.2- Obviously these functions verify the Helmholtz equation (1) as a hnear 
combination of solution, and with the same eigenvalues. Finally, one has to 
note that since T is an invcrtible matrix, only 0„ = would give ipn — 0. 
This proves the isospectrality between M.i and A42- 



2.3 Trajectories on the translation surfaces 

The matrix formalism introduced in subsection 2.1 allows us to construct a 
"movement matrix" M from any path drawn on the translation surfaces Aii 
or M.2, adapting the method explained in [13]. Each tile has 6 neighbours, 
and any path is drawn on a sequence {ii,i2, ■■■,in) of tiles such that ik+i is 
a neighbour of ik, provided it does not hit any vertex (we call "vertex" a 
point on the surface where 4 tiles join, or a point in the middle of a hor- 
izontal edge; there arc 4 scattering vertices, the others are non-scattering 
ones). Let us call Ak G {R,L, U , U , D , D} the matrix corresponding to 
the movement from ik to ik+i, i-e. the matrix verifying (^fc)ife,ife+i = 1- If 
we define the movement matrix M as nfc=i^fc' ^^^^ verify Mi^^^i^ — 1. 
As a product of permutation matrices, M is a permutation matrix, and 
it maps tile ii onto tile in- For instance the path drawn in Figure 3 cor- 
responds to a sequence of tiles (1,6,4,4,6,5) and to a sequence of ma- 
trices M — RU LU R. Reciprocally, for any product M — AiA2...An of 
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Figure 3: A path drawn on A^i 
matrices belonging to {R,L, U ,U , D , D}, if Mjj = 1 then there is a se- 
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quence {ii = i,i2, ■■■,in = j) of tiles such that ik+i is a neighbour of ik and 
{Ak)ik,ik+i ~ Note that infinitely many sequences (ii,i2, ■■■,in) give rise to 
the same movement matrix, because there is an infinite number of sequences 
and only a finite number (7!=5040) of permutations. 

Let us now consider a sequence of matrices {A''^\ A^2 \ •••) An^), and the move- 
ment matrix M^^'^ = A^-^^ A2^\..An^ ; then for any i it defines a unique sequence 
of tiles (i = il, ^2, in) such that (^fc^^)ifc,ij;+i = 1- This sequence has the 
property that M^^^-^ — 1. According to the relation (7), each ^4^^^ verifies 
^{i) ^ 7--i^(^2)^^ ^^^^^ implies that the matrix M^^) = A^^^ Af\..A^n^ veri- 
fies 

TM^^) = M^^^T. (10) 

The interpretation of this commutation relation is the same as the interpre- 
tation of (7) for the individual A^ (see Equations (8) and (9) with R replaced 
by M): if one can go from ii to i„ by a sequence of tiles glued in a specific 
way, then the three pre-images of i„ are the tiles obtained by starting from 
the three pre-images of ii and following a sequence of tiles glued in exactly 
the same way. To any path drawn on not hitting corners it is possible to 
associate the sequence of tiles on which it is drawn, starting from a tile i and 
finishing on a tile j; it is therefore possible to draw an identical path starting 
from any of the three pre-images of i; this path will necessarily arrive to one 
of the three pre-images of j. 

2.4 Equality of the length spectrum 

The two surfaces A4i and M2 are pseudo-integrable, therefore their periodic 
orbits occur in families of parallel orbits of same length [16]. It has been 
shown in [13] that there is an exact one-to-one correspondence between the 
pencils of periodic orbits of A^i and those of Ai2- The main argument is 
the principle explained in section 2.3 that for any path drawn on one of 
the surfaces there is a corresponding sequence of tiles and a corresponding 
sequence of matrices. The product of these matrices, the "movement matrix" 
M, has the property that Mij = 1 if and only if the sequence of tiles goes 
from i to j. For any periodic orbit on the last tile has to be equal to 
the first one, and a closed path going from tile i to itself has a movement 
matrix verifying M^^^ = 1. The quantity tr(M*^^)) is therefore the number of 
tiles from which one can start and come back to oneself after a sequence of 
tiles giving the movement matrix M^^\ Since the commutation relation (10) 
implies that tr(A'/*^-'^^) = tr(M*^^)), there is the same number of tiles having 
this property in 7^2- So for any periodic pencil drawn on this sequence 
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of tiles there will be tiM^^^ identical copies of it (in particular with same 
length and same width) on A^i and the same number on A^2, hence the 
bijection between the periodic orbits. As an illustration, Figure 4 shows the 
two pencils of periodic orbits in a given direction on A4i and A42: the grey 
orbit has same length and same width on both surfaces, and so does the 
white orbit. 




3 The Green function 

Relation (3) implies that the imaginary part of the trace of the (retarded) 
Green function is identical for two isospectral billiards. Here we will be more 
precise and express the Green function of AI2 in terms of the Green function 
oiMi. 

3.1 Tile modes and normalization 

Each surface M.i and A^2 is made out of 7 tiles glued together. We will 
call tile modes the solutions of the Helmholtz equation (1) on a tile with 
periodic boundary conditions, that is a torus of size u x v. The eigenvalues 
corresponding to these tile modes will be denoted E^, and the corresponding 
eigenfunctions Xt- 
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There is a one-to-one correspondence between these tile modes Et and 
a subset of the spectrum £"„ common to the surfaces A4i and A42- First, 
any eigenfunction xt corresponds to a solution ipt{a,i) (or (f)t{a,i)) of the 
Helmholtz equation on A^i (or M.2) by simply taking the function equal 
to (or proportional to) Xt{o) on each tile: the periodic boundary conditions 
for the tile eigenfunctions Xt will make tf^t and 4>t continuous. In order to 
normalize correctly to 1 these eigenfunctions, one has to set 

V't(a, i) = 0t(a, i) = -^Xt{a)- (H) 

Reciprocally, for any eigenstate 0„ of the surface A^i, the function defined 
on a tile by 

7 

Xn{o) ^^(t>n{a,i) (12) 

i=l 

is either the function or an eigenstate Xt of the tile. We will use the 
subscript s when the state 0„ verifies the condition 



(a,i) = (13) 



and the subscript t when it docs not. The set of cigcnstatcs of M.i (and, in 
the same way, the set of eigenstates of M.2) can then be partitioned into two 
sets: the eigenstates 0^ which are also the eigenstates of the tile and do not 
have the property (13), and the eigenstates 0s which have the property (13) 
(pure "surface" states). 

The normalization constant An in (5) will take a different value whether 
the function 0„ belongs to the set of {(pt} or {0s}, as we will see now. Equa- 
tion (5) expresses each eigenfunction ■j/'n of as a sum of an eigenfunction 
0„ of All with the same eigenvalue, taken at three different points of M.i 
(or rather at a similar point on three different tiles). The normalization 
condition on ip can be written 



1 



\il)n{,x)\^ dx = ^ \ipri{a, k)f da 

M2 f^—i "'tile 

7^ 5Z Tk,iTk,j / 0n(a, i)(l>n{a, j)da 

•^n j j^i ^tile 

= -H^^lj I 'Pn{(^^i)'Pn{a,j)da (14) 
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where we have used the fact that T is symmetric and summed over k. Since 
the matrix is equal to {T'^)i,j = 1 + 25ij {5ij is the Kronecker symbol), we 
get 

= XI / 4>n{a,i)(l)nia,j)da + 2'^ / \(f)n{a,i)\^da 

; .-1 JtWe .-I ^tile 



tile 



i=l 



da + 2 |0n(fl, i)\ da. 

'Ml 



(15) 



For the tile modes ipt, we can replace (p by x/V7 from Equation (11). This 
yields 

At = — \xt{ci)fda + — \xt{a)fda = 9 (16) 

' Jtilc ' Jtilc 

because the functions Xt ai'c normalized to 1 on the tile. For the non-tile 
modes (ps, which verify Equation (13), Equation (15) gives — 2. Finally 
we have 

V^t (a, ^ XI '^hj^riia, j) (tile modes) 

3 

ips{(i,i) = —^^Tij(f)n{a,j) (surface modes). 

j 

(17) 



3.2 The Green function of the translation surfaces 

In this section, we want to express the Green function on the surface in 
terms of the Green function of the surface M.i. We will use the expansion 
over eigenstates for the advanced Green function of 



E - En + ie 



(18) 



where ipn and En arc respectively the eigenfunctions and the eigenvalues of 
/A2- We have to split the sum over n into a sum over the tile modes ipt and 
the non-tile modes ips, and replace ip by its expression (17): 



G^'{a,i:h,i) = -y 



1 T.i',j'Ti,i'Tjjf(f)s{a,i')(j)s{b,j') 



E-E, + ie 



(19) 
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If we add and subtract (1/2) we get 



-- E / E '^^^"'^^^'^^^^''^ (20) 

i'j' t * 

The first sum is a sum over Green functions of A^i; in the sum over the 
remaining t modes, we replace (pt by its value given by Equation (11). This 
gives 

G^'\a,i;b,j) ^]^Yl^,,>T,,,G^^\a,i'-h,3') - 1 J] r,,,r,-,,GW(a; ft) (21) 

where G^^\a;b) is the Green function on the tile. This Green function does 
not depend on i' or f: therefore the sum over i' and / can be performed: 
Yli' j' — 9 since each row and each column of T has three I's. Finally 

the Green function of the surface A^2 is 

G^'\a,t-Aj) = lY.T^,^'TJ,fG^'\a,^'■Af) - V)(a;6). (22) 



Obviously, a similar relation can be obtained that expresses as a sum 
over functions G^'^\ 



4 Diffractive orbits 

4.1 Stovicek's formalism 

We know from [11] that the exact expression for the Green function can be 
obtained for billiards with a polygonal enclosure as a scattering series, where 
each term is a sum over scattering paths made of straight lines and scatters 
on the singular corners. Each scattering path contributing to G{a, b) is made 
of a starting leg of length ro going from the initial point a to a scatterer, 
then an alternating series of diffractions with angles G M followed by legs 
of length Ti, 1 < i < n, the last leg of length r„ going from the last scatterer 
to the final point b. The expression given in [11] for scattering on straight 
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reflectors can be adapted here and gives 



^ ' n vertex 



— OO 



'OO 



dsids2---dSnHQ^^ [kR{si, S2, 



paths 



n 



X n 



— 71- 



■2 ' 



(23) 



k=l 



where 



R^si, S2, Sn) = (ro + rie^i + r2e''+'' + ■■■ + r„e''-^+^^+-+^") 

X (ro + rie-'i + rse-'^-"^ + ■ ■ ■ + rne"'!-"^— . (24) 



The A;-th diffraction angle (fk is equal to Mk'jk + ^^fc; here 7^ is the measure 
of the angle at the singularity {An or Gvr in our case), Mk is the winding 
number (i.e. the number of times the path winds around the singularity), 
and < ^fe < 7fe. A schematic example of such a scattering orbit is provided 
at Figure 6. 

4.2 Diffractive orbits 

Let us now look at the "saddle-connexions" [17] which we define as the 
geodesies joining two diffracting vertices. If u and v are the width and the 
height of a tile, the vertices (scattering and non-scattering) are in the direc- 
tions {mu/2,nv) with m, n G Z. Each pair {171,11) with m and 11 co-prime 
will give the direction of a saddle-connexion. If we define 



for p = {m, n) E Z"^, m and n co-prime, the lengths of the saddle-connexions 
will be of the form kip, k e N. These saddle-connexions will be the legs of 
the scattering paths in the expansion (23). 

But there is no one-to-one correspondence between the saddle- connexions 
of A^i and those of M.2- In fact, there are saddle- connexions in M.2 that 
cannot be found anywhere in Ali. For instance, the dashed saddle-connexion 
on the surface AI2 in Figure 4 in the direction (m = 1, n = 2) has a length 
4lp, which is twice as long as any of the saddle- connexions existing in A^i. 
A way of understanding how this difference between the sets of diffractive 
orbits still allows isospectrality is to prove the relation (22) between the Green 
functions, using their expression (23) as a definition. This will be done in 
the following section, by establishing a certain correspondence between the 
scattering trajectories. 



(25) 
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4.3 Correspondence between difFractive orbits 

Let us consider a contribution to (23) for any of the two surfaces. It is a suc- 
cession of scatters and saddle- connexions of lengths {kpjp^, kpjp^, kpjp^), 
where the ki are integers and the k are defined by (25). In the case of forward 
diffraction (9^ — n [27r]), the integral over Sk has a pole at Sk — 0. But the 
discontinuity arising from this singularity is canceled by an opposite discon- 
tinuity in the term of order (n — 1) in (23), as it should for physical reasons 
of continuity. This cancellation is discussed in [11]. This means that in the 
expression (23) one has to interpret any term containing a forward scattering 
9k — TT [27r] between a and b as the limit for e — > of a term corresponding 
to a path from to b^, with — > a and b^ b. This shifted path from 
to feg will have two contributions: a straight path from to b^ missing the 
singularity, plus a sum over all the scattering contributions starting from 
and winding any number of times around the singularity (see Figure 5). 




Figure 5: A contribution to the Green function in case of forward diffraction 

Furthermore, this still holds even if there is no scatterer at the vertex 
because in that case the series of diffractive terms adds up to zero: when 

7k = 27r, 

°° 27r 

E (^,Mfe + 7r + isfe)2-7r2 ^^^^ 

Therefore any saddle-connexion of length kip can be replaced (as in Figure 
5 but imagining now that a and b are scattering vertices and that the x 
corresponds to a non- diffracting vertex), by a sum over all possible paths 
made of straight lines of length A';,;/p parallel to the saddle-connexion and by 
windings around the non-scattering vertices any number of times, with the 
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condition that '^ki = k. So, any contribution to (23) can be decomposed into 
an infinite sum of paths made of a fixed number of windings around scattering 
vertices and any number of windings around non-scattering vertices. This is 
illustrated at Figure 6. 



b 




Figure 6: A contribution to the Green function in case of forward diffraction. 
The filled circles are scattering vertices, the empty ones are non-scattering 
vertices. 

An orbit on the torus going from a to 6 consists only of a straight path 
between a and b, since there is only one vertex and it is non- diffract ive. But 
again, one can decompose this straight path into a sum over paths going 
from the vertex to itself and winding any number of times around the vertex. 
These contributions will add up to zero. 

Let us now fix two points a and b on the torus and consider a path of 
the type described at the left-hand side of Figure 6, with any given num- 
ber of windings around the vertices. All the contributions to G^'^\a,i]h, j), 
G''^\a,i;b, j) or G^^\a;b) are of this type. Reciprocally, one can write each 
Green function as a sum over such contributions weighted by 1 if such a path 
exists between {a,i) and {b,j) and by otherwise. 

Since such a path does not hit any vertex, the sequence of matrices 
(^1, ...,An) corresponding to the movement is well-defined (see section 2.3), 
and one can construct the corresponding movement matrix M = AxA2...An- 
This path always exists on the torus, which means that it is always a contri- 
bution to G*^*^(a, b). According to section 2.3, it will exist on the surface M.v^ 
1/ = 1, 2, if and only if a and h are on tiles i and j such that M^- = 1. There- 
fore any such path contributes to the right-hand side of Equation (22) with 
a weight -1/2 coming from GW(a,6) and with a weight Y,i'^/'^)iTi,i'Tj,j'), 
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where the sum runs over all the pairs i', j' such that the orbit exists between 
tile i' and tile j', i.e. such that M-pj, ~ 1. The total weight associated to 
this path in the right-hand side of Equation (22) is therefore 

\T.T^,i'Tj,fM^>-\. (27) 

Using the commutation relation (10), we have 

T.T^,^'TJ,J'M^Py=iM^'^T')... (28) 

But it can be computed directly from the expression (4) for T that {T'^)ij — 
1 -|- 2Sij. Therefore 

(^^'^n,, = E<?+2<). (29) 

k 

Since M^^^ is a permutation matrix, the sum over a column is equal to 1. We 
get, from equations (28) and (29), 

E^.^'^.^'^y = l + 2<?- (30) 

The weight (27) of a path in the right-hand side of Equation (22) is therefore 
equal to M^'^j , which is exactly its weight in the expression of G^'^\a,i] b,j). 
This is an alternative way of proving the equality (22), using the formula 
(23) for the Green function instead of the isospectrality of the surfaces. 

In fact. Equation (30) provides a more precise result: it states that 
there is a "one-to-three" correspondence between the diffracting orbits of 
Ai2 and those of A^i. The main point is that for given tiles on Ai2, 

ji'^i^i'Tj. j'M^}^-, is the number of times a path of movement matrix M^^^ 
appears among the 9 Green functions G^^^a,i';b, j') with i' pre-image of i 
and f pre-image of j. So if a path exists on AI2 between tile i and tile 
j (which means Alj^'^j = 1), then 3 copies of it exist in A^i (and this will 
necessarily be between the 3 pre-images of i and the three pre-images of j), 
whereas if a path does not exist in Ai2, then by Equation (30) only 1 copy 
of it exists in A^i (between one of the three pre-images of i and one of the 
three pre-images oi j). 
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5 Back to periodic orbits 



A final question one might want to ask is the following: why is there a one- 
to-one correspondence between periodic orbits, whilst each periodic orbit in 
Ai2 should correspond to 3 orbits in Aiil The answer is that these 3 orbits 
are not necessarily periodic. Consider a periodic orbit in Ai2, not hitting any 
vertex (almost all of them verify this condition), going from tile i to itself 
with a movement matrix M^'^\ Let us call {^1,^2,^3} the three pre-images 
of i. Then, according to what we said in the previous section, there are 3 
copies of this orbit in Aii going from k G {^1,^2? ^^3} to k' G {^i,^2)'^3}- But 
these copies in Aii are not periodic if k 7^ k'. Moreover, we know that there 
are in fact 7 copies of this orbit in Mi, if we do not restrict ourselves to 
the pre-images of i and j as starting and ending tiles. Therefore for periodic 
orbits the correspondence is more global: the condition of periodicity imposes 
that we take into account not only orbits from pre-images to pre-images, but 
orbits from any tile to any tile. Then, as we already said, the number of such 
orbits is trM^-'^^ which is equal to trM^^^ We should therefore speak of a 
trM^-'^^-to-trM^^^ correspondence between periodic orbits, the value of trM^^^ 
depending on the periodic orbit considered. 



6 Conclusion 

For the pair of billiards given in Figure 2, which are made of rectangular 
tiles glued together, there exists a set of neighbour matrices {R^^\ L^^\ 

describing the way these tiles are glued together. 
The essential feature which accounts both for isospectrality and for the corre- 
spondence between paths on the surface is the existence of a "transplantation 
matrix" T which has 3 properties: 

• it is invertible (otherwise one of the spectra would just be a subset of 
the other) 

• it is not a permutation matrix itself (otherwise the two domains would 
just be congruent) 

• it has the commutation property TM^^^ — M^'^^T for all neighbour 
matrices M, which assures that smoothness at all the segments between 
the tiles and boundary conditions will be satisfied. 

It is therefore possible to generahse the previous analysis to all the pairs 
of isospectral biUiards made out of base tiles glued together, provided one can 
find a transplantation matrix T having these 3 properties. It turns out that 
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all the examples of pairs of isospectral billiards, as far as is known (see [6] 
and [18]), are constructed by the same application of a theorem by Sunada 
[19] and consist of tiles of a given base shape glued together. Moreover, 
Sunada's theorem implies the existence of a transplantation matrix in each 
case [18]. Therefore Equation (28) always holds, and the same arguments 
can be adapted to establish a correspondence between diffractive orbits in all 
these other known cases. However, the correspondence depends on the entries 
of the matrix T and T^, and might be more comphcated in the general case. 
The same results could also be worked out for pairs of isospectral billiards 
which would not be based on Sunada's theorem but would nevertheless have 
a transplantation matrix. It is not known however if such billiards exist [18]. 
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A Appendix 

The matrices that describe the gluings between the plates of the translation 
surfaces M.\ and are obtained by Figure (2). For the structure M.i they 
read 
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For the structure M.2 they read 
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(32) 



One can verify exphcitly that the commutation relations (7) hold for these 
matrices. 

It is of some mathematical interest to notice that the groups T^^^ and F^^) 
generated by the neighbour matrices L(^), C/M, C/M, ^H, ^M} for 

1/ — 1,2, are subgroups of 5*7 (the group of permutations of 7 elements) of 
order 168. These groups turn out to be isomorphic to the linear group ^2(7) 
(also called PSL{2,7)), which is the group of automorphisms of the finite 
projective plane of order 2, or Fano plane (see [20] for a definition). The 
matrix T is nothing but the incidence matrix of the graph corresponding to 
the Fano plane. 
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